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Abstract 

Wc define the notion of a spanning tree generating function (STGF) ^ a„z", wliich 
gives tfie spanning tree constant when evaluated at z = 1, and gives the lattice Green 
function (LGF) when differentiated. By making use of known results for logarithmic 
Mahler measures of certain Laurent polynomials, and proving new results, we express 
the STGFs as hypergeometric functions for all regular two and three dimensional 
lattices (and one higher-dimensional lattice). This gives closed form expressions for 
the spanning tree constants for all such lattices, which were previously largely un- 
known in all but one three-dimensional case. We show for all lattices that these can 
also be represented as Dirichlet L-series. Making the connection between spanning 
tree generating functions and lattice Green functions produces integral identities and 
hypergeometric connections, some of which appear to be new. 



1 Introduction 

There is a simple connection between spanning trees (ST) and lattice Green functions 
(LGFs). This connection has been used previously, as detailed below, but has not pre- 
viously been systematically exploited. It may be used to derive the ST growth constant 
from the LGF, and as a byproduct, we derive certain integral identities involving com- 
plete elliptic integrals of the first kind, and connections between different hypergeometric 
functions, some of which appear to be new. 

^email: tonyg@ms.unimelb.edu.au 
email: mathewrogers(§gmail . com 
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In the next section we define spanning trees and spanning tree generating functions 
(STGFs) and briefly review their properties, particularly those relevant to this study. In 
the following section we define and provide relevant reviews of properties of lattice Green 
functions. In the fourth section we define logarithmic Mahler measures, and give some 
results needed in subsequent derivations. In section 5 we exploit the connection between 
LGFs and STGFs. In particular, the integrals defining the STGFs are readily seen to 
be expressible in terms of logarithmic Mahler measures. From the existing literature on 
Mahler measures [32,33], plus significant extensions given here for the simple-cubic case, 
all of the STGFs for the standard two- and three-dimensional lattices are shown to be 
expressible in terms of hypergeometric functions^. We also express these constants in 
terms of Dirichlet L-series. In Section 6 we use these results to produce integral and 
hypergeometric identities. Section 7 briefly outlines the connection between spanning 
trees, dimer coverings and the Ising model. 

2 Spanning trees on a lattice 

A spanning tree on a graph G is a loop-free connected graph connecting all sites of G. 
The number of spanning trees on a graph G we denote Tq- For a class of graphs called 
recursive^, which includes regular lattices, Shrock and Wu [37] have proved that the num- 
ber of spanning trees grows exponentially with the number of sites of the lattice. That 
is to say, for a regular lattice C oi N sites, the number of spanning trees on C, denoted 
Tc{N) grows like e"^^ for N large. The limit limAr_5.oo jf ^ogTc{N) = Xc exists, is greater 
than zero, and depends on the lattice C. We call this limit the spanning tree constant. 
Lyons [29] refers to it as the tree entropy. 

The exponential growth with the number of sites is not the case for all graphs. For 
example, as pointed out in [37], for the linear chain of sites there is but one spanning 
tree, whereas for the complete graph of N sites the number of spanning trees grows faster 
than exponentially. 

There are several ways of calculating the number of spanning trees for a given graph. 
The standard graph-theoretical method is by construction of the Laplacian matrix, see, for 
example [5]. Spanning trees can also be related to a special value of the Tutte polynomial. 

As pointed out by Wu [45], if the underlying graph is periodic, that is to say a lat- 
tice, then following Fortuin and Kasteleyn [14] Tc{N) can be expressed as the partition 
function of a lattice model. This representation then allows one to calculate Tc(N) as the 
partition function of an ice-type model on a related lattice [1]. In this correspondence, the 
partition function is evaluated as a Pfaffian, and the spanning tree constant is given as a 
c/-dimensional integral. This integral may be written as 



^ After the completion of this work we became aware of the recent paper of G S Joyce [23], who studied 
the problem on the face-centred cubic lattice and obtained results in agreement with ours, as well as several 
additional results on related problems 

^ A class of graphs is recursive if it can be constructed by sequential addition of a given subgraph. Thus 
regular lattices are recursive. 
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where q is now the co-ordination number of the lattice^ and A(£) is the structure function of 
the lattice. The co-ordination number is just the number of nearest neighbours per lattice 
site, while the structure function is the Fourier transform of the discrete step probability 
function. For example, for the d-dimensional hypercubic lattice C^, one has q = 2d and 

A(Cd) = ^(cos ki + cos k2 + ■ ■ ■ + cos kd). 
d 

The above integral expression for the ST constant has subsequently been independently 
derived by graph theorists, first by Burton and Pemantle [8], and in a slightly simpler 
form by Lyons [29]. 

There is a considerable body of literature on the problem of evaluating for a variety 
of lattices C In two dimensions this can usually be done exactly, but not generally in 
higher dimension (the hyper-body-centred cubic lattice is an exception to this remark). 

In two dimensions the results for the square, triangular, honeycomb were first given 
by Wu in [45] and by Shrock and Wu [37] for the kagome lattice. They are: 

^sq = / / d/ci(ifc2log[4 - 2(cos/ci + cosA;2)] (1) 

4 111 AG 
= -(l-^ + ^- ^ + ---) = — = 1.166243616 . . . 

TT 3^ b'^ 7^ TT 



1 

J J d/ci(iA;2 log[6 — 2(cos fci + COS + cos(A;i + ^2))] 
^(1-^ + ^-^ + ^ + ...) = 1.61532973... 



^honey ^tri/2 
^kagome = ( Af ri + log 6) /3. 



We remark that these results can also be expressed in terms of primitive Dirichlet 
L-series [44]. In particular, 

Xsq = -L^i{2) (2) 
vr 

Atri = ^^ -^-3(2) 

A - 15^/3 

Hagome = | ^^L_3 (2) + ^ log 6 j . 



In the above equation, G is Catalan's constant; the result for the honeycomb lattice is 
a consequence of a theorem [37] linking \c on a given lattice C to A^* on the dual lattice 
C*. A variety of other two-dimensional lattices have also been considered in [37] (as well 
as some lattices of higher dimension), while in [11] it is pointed out that there is more 

^Not to be confused with the number of states of the Potts model 
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than one way to define a lattice, wliich gives rise to a range of different integral expressions 
for the spanning tree constant. In this way Chang and Wang found a number of integral 
identities, by virtue of different choices of integrand, which all evaluate to the spanning 
tree constant for the lattice in question. In [10] some non-regular two-dimensional lattices 
are considered, and a systematic treatment of such latices is given, including the evaluation 
- in most cases numerical rather than exact - of the spanning tree constant for a wide 
variety of such lattices. 

In what follows, it will turn out to be useful to generalise the integral and form what 
we call a spanning tree generating function (STGF). We define it as 

Tc{z) ^ log q+j^ j\h--- dkd\og[l/z-K{k)]. (3) 

Clearly, the spanning tree constant A£ = Tc{^)- It immediately follows that 

- - (2^ £ • • • £ *rr^ - ^^<»- 

The integral in eqn(4) will be recognised as the lattice Green function of lattice C. 
For some lattices, such as the honeycomb and diamond lattices, the LGF is given by 

^^<»'^> = P^£'*-£*^TriW' 

In that case the STGF is 

Tc{z)=\ogq+]^-^^ J\h--- J\ka\og[l/z^-Aik)], (6) 

and eqn (4) still holds. 

This generalisation is not new. It appears to have first been used by Rosengren [35], 
who obtained an expression for the spanning tree constant on the simple cubic lattice by 
making use of known recurrences for the LGF of the simple cubic lattice, and integrating 
these. More recently, Joyce, Delves and Zucker [24], inspired by an integral occurring in 
work of Baxter and Bazhanov (unpublished), which in structure is of the form of eqn(3), 
instead investigated the derivative function, which gives rise to an integral of LGF type, 
as in eqn(4). In [22] Joyce studied integrals equivalent to the STGFs, with a view to 
obtaining accurate estimates of the spanning tree constants. He evaluted the STGF for 
the body-centred cubic lattice in terms of hypergeometric functions. He also evaluated the 
spanning tree constants of the three standard three-dimensional lattices to extraordinary 
accuracy, nearly 200 decimal digits. Further, in 2005, Glasser and Lamb [15] considered 
the lattice spanning tree entropy for various two-dimensional lattices, and used it to derive 
some anisotropic triangular lattice Green functions by differentiation. 

However here we exploit this connection in a systematic manner. The existence of 
known exact results for some lattice Green functions can now be integrated and provide 
new, simpler, representations for the spanning tree generating function, as well as new 
integral identities. More precisely, we have 

Tc{z)= log q- [ ^^dz. (7) 
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Note first that we have restored the constant of integration, logq, lost in the differ- 
entiation. Secondly, when the indefinite integral can be evaluated, and the STGF found, 
the spanning tree constant can be found by evaluating the STGF at z = 1. In most cases, 
particularly for lattices of dimension greater than 2, the integral cannot be performed. In 
that case we can still get an expression for the spanning tree constant. First, note that, by 
definition, the LGF can be written -P(0, z) = 1 + 0{z). The integrand of (7) therefore has 
a simple pole at the origin, which contributes a term log z to the STGF. This of course 
vanishes when evaluating the STGF at 2; = 1, while all higher order terms vanish at z = 0. 
So if we artificially remove this pole at the origin of the integrand and integrate 
between and 1, we will obtain the spanning tree constant. That is to say, 

Xc = logq- dz. (8) 

Jo z 

This result has been previously obtained by Lyons [29]. 

In the next section we review some results for lattice Green functions which will be 
needed in our subsequent development of spanning tree generating functions. 



3 Lattice Green functions 

For a translationally invariant walk on a d-dimensional periodic Bravais lattice, a natural 
question to ask is the probability that a walker starting at the origin will be at position / 
after n steps. The probability generating function is known as the Lattice Green Function 
(LGF). It is 

p(i-z)- ^ r r^MziMl 

^(^'^)-(2^y„. i-,A(fe) ■ 

So that [z^]P{l; z) is the probability that a walker starting at the origin will be at I after 
n steps. A(/c) is the structure function of the lattice walk, as noted above. 
The probability of return to the origin is 1 — 1/P{0; 1). 

, 1 r r dkidk2 
P(0;1) ' ' ^ ^ 



(2^)2 1 _ A(fc) 

Since P{0; 1) diverges for two-dimensional lattices, this leads to the well-known result that 
the probability of return to the origin in two dimensions is certain. In three dimensions 
these are the celebrated Watson integrals. The history of their development and evaluation 
has recently been authoritatively given by Zucker in [47]. 
Of broader interest are LGFs defined by: 

1 r r d'^k , , 
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For the regular two-dimensional lattices the structure functions are^ 

Hk)honeycomb =^ (1 + 4 COS^ X + 4 COS X COS y) , 
M^)square =^(cOs/ci + COsA;2), 
Mk)triang =^(cOS ki + COS k2 + COS {ki + k2)). 

The corresponding LGFs are: 



P{^;z)honev = =K(fc) (11) 



where 



(3-z)V^(3-z)(l + 2)' 

where K(2;) is the complete elliptic integral of the first kind. For the square lattice, the 
result is remarkably simple, 

P(0;z),, = -K(z), (12) 

TT 

while for the triangular lattice the LGF is: 

P(0;z)t„ = ^K(A;') (13) 

where c = (a + 1)(6 — 1), and 



3 / 6 3 / 6 

a = - + 1 - Js + -, and 6 = - + 1 + Js + - 

Z V 2 V -Z 



and 



^,_,/2(6-a) 



c 

For the square lattice, we can also use the equivalent structure function 

^{k)square = COS ki COS k2. 

as the square lattice can be considered as the two dimensional hyper-cubic lattice, which 
gives the first form of the structure function (above), or as the two-dimensional hyper- 
body-centred cubic lattice, giving the second form. While the integrands are clearly dif- 
ferent, the integrals are equal, that is to say, 

^ ^ r r dkidk2 ^ _j_ r r dkidk2 

[U, Z)square ^^^y ^ _ , ^^^^ ^ ^^^y ^ _ ^^^^^ • 

(14) 



^Because the honeycomb lattice has two types of site, the expansion parameter z in eqn (10) should be 
replaced by z^. 
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Similarly, for the honeycomb lattice we can exploit the duality with the triangular lattice 
and write the structure function as 



2 1^ 



It follows that 

-P(0, Z^honey 



1 



(27r)2 
1 



dkidk2 



-77 J —71 1 
71 /"TT 



[I + 3 [cos ki + cos k2 + cos(A;i + ^2)]) 
dkidk2 



(15) 



(2vr)2 J^j^J^T, 1 - + 4cos2A;i +4cos/ci cosA:2)' 



The LGFs for the three-dimensional lattices are also known. For the simple cubic 
lattice one has: 



PiO;z) 



1 



{tt) 



JO 



dki dk2 dk^ 



1 — f (cos ki + cos k2 + cos ks 



Joyce [21] showed that this could be expressed as 



1-9^^ 



where 



kl 



(1-0^(1 + 30 
16^3 



-Kf 



TT 



(1-6^(1 + 30^ 



with 



e = (1 + Vl- ^2)^1/2(1 _ ^l_^2/9)l/2_ 

For the body-centred cubic lattice one has: 



P{0-z) 



(7r)2 



JO 



d/ci dk2 dk^ 



1 — z(cos ki cos A;2 cos /ca 



Maradudin et al. [30] showed that this could be expressed as 

1 2 



P(0;z) 



-K(A:2) 
vr 



where 



A:2 = i-iyr3^ 

^ 2 2 ^ 



For the face-centred cubic lattice one has: 
P(0;z) ^ 



Jo 



TT /'TT /"TT 



JO 



dki dk2 dk^ 



1 - |(ciC2 + C1C3 + C2C3) 



where Cj = cos/cj. Joyce [21]) showed that this could be expressed as 



2\2 



P{0;z) 



(1-C)3(1 + 3C) 



vr 
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where 



and 

e = (1 + vT^rH-i + vTTsJ). 

Finally, for the diamond lattice one has: 

1 r r r dkidk2dk3 



P{0;z) 



i'^)^ Jo Jo Jo 1- ^(l + CiC2 + CiC3 + C2C3)' 

which Joyce [20] pointed out gives 

P{6;z) = ^K{k+)K{k.); 

where 

kl = l±lz\A-z^f/^^-\{2-z^){l-zT^'\ 

The derivation of these results and some history and erudite discussion can be conve- 
niently found in the book by Hughes [19]. 

4 Mahler measures 

The (logarithmic) Mahler measure of an n- variable Laurent polynomial is usually defined 
by 



m(P(2i, . . .,Zn)) :-- 



C ■■■ C log |P(e2^'^i, . . . , e2^'^")|(i0i ...dOn. 
Jo Jo 



In the two-variable case, we will be interested in the following polynomials, which we 
denote as shown: 

m{k) := m ( k + X -\ h y H — ), (16) 

\ X yj 

n{k) := m{x^ + + 1 — kxy) , 

g{k) := m{{l + x){l + y){x + y) - kxy). 



The Mahler measures above can be represented in terms of hyper geometric functions. 
The generalized hypergeometric function is defined by 



{ai)n . . . {ap)n 
(&l)n.-.(&g)n uV 
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where {a)n = r(a + n)/T{a). It was proved in [31] and [28] that 

m{k) := JR(log(A:)-^4F3(%'^^i|)), (17) 

2 / 1,1,1,1:27 



n{k) := ^llog{k)--.^F- 



3'3' 



y^S'* I 2,2,2 f^3 



/(4 + fc)(fe-2)n 2k^ ^(l±ii.27k^ 



8k „ /s 4 1 , 27A; 



3'3 



,1,1; 



{k-2)^''^^ V ^'^'^ (fc-2)3 



The derivation of the formula for g{k) requires a modular expansion due to Stienstra [39], 
as mentioned in [28]. 

In the three-variable case we will be interested in the following polynomials: 



p{k) := m(yfc+(x + i) (y + i) (z + i)), (18) 

^ 1 1 1 

s[k) := m I k + X A |-yH \- z ^ — 

V X y z 

m := ni(4-^+(x+i)(,+ i) + (. + l)(z + l) + (, + i)(z+l 

In [32] one of us proved that for u sufficiently large 

P(^) := \f2{u), (19) 

^ 1 , /9(3 + nY\ 3 ^ 9(3 + n-Y 
s(3(n + u )) := — jaI t. H /4 s 

1 /(16-n)3\ ^ 8 / (4-^x)3 



where, for |n| sufficiently large, 

h{u) = K(Mu)-V4(ti4'y^^)), (20) 

U{u) = Kflog(n)-^5F4te''^^''' 



u 



2,2,2,2 



For our purposes here, we point out that 
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s{l/u) 



log(u) 
log(u) 



2n /2j^ 

n 



n=l 
n=l fc=0 



fc=0 



2k 
k 



2n - 2k 
n — k 



2k 

k 



for \u\ < -, 
6 



for lul < — . 
' ' 16 



(21) 



5 Connections between spanning tree generating functions, 
Mahler measures and lattice Green functions. 

From equation (4), one can clearly calculate the spanning tree generating function from 
the lattice Green function (9), up to a known additive constant, notably logg, where q 
is the co-ordination number of the lattice. Furthermore, from the results given above for 
Mahler measures, we can in principle express the STGFs for all lattices in terms of Mahler 
measures, and hence in terms hypergeometric functions. 

5.1 Two-dimensional lattices 

For example, for the square lattice, from eqns (16) and (17) one has 



log 4 + 



1 

4^ 



log 



1 1 



(cos ki + cos k2) 



log4 + ^ / / log 
m(4/z) 



dki ■ dk2 
dki ■ dk2 



log 4 - log z - ^4i^3 (^'2 2^ ; 



(22) 



The second equality follows since the function inside the logarithm is greater than or equal 
to zero for < z < 1, allowing us to introduce the absolute value sign inside the logarithm, 
thus turning the expression into a logarithmic Mahler measure. If we substitute x = e^'^^^^ 
and y = e'^'"^^^, then the third equality immediately follows. 

In this case, we can also obtain the result by direct integration of the LGF, since for 
this lattice P{0;z) = ;|K(z), so from eqn (7) 

T,,{z) = log4 - ^ 1 ^dz = log 4 - log . - ^4^3 (^isV ; • (23) 
From [33] the r.h.s can be further simplified to 
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It immediately follows that, 



^sq 



log 4 



1 2 



^K(z) - 1 



z 



dz = log 4 — -^4^3 



3 3 
2 ' 2 • 1 
2,2,2 ' 



1,1 



(24) 



3-f2 1 3 , 



4G 
vr 



For the hexagonal lattice, 



Thex{z) =log3 + ^ 




log 



— 7T J —7T 



-TT (l + 4cos^ ki + A cos ki cos ^2) 

9 



87r2 



log 



Let us set 



-2 — 1 ) — (4 cos^ fci + 4 cos k\ cos ^2) 
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dA:i(i/c2 (25) 
dkidk2 (26) 



= ^ - 1. 



The function inside the logarithm is greater than or equal to zero if A; > 8, or equivalently 
if —1 < z < 1. This allows us to introduce an absolute value sign inside the logarithm. 
We obtain 

1 

Thex{z) =-^-2 J J log j/c — (4cos^ /ci + 4cos A:i cos ^2) I 'iA:i(iA;2 

=- / / log |A; — 4cos (27rA;i) (cos(27rA;i) + cos (27rA;2))| (iA;idA;2 
2 Jo Jo 

= }- I I log\k — 8cos {Airki) cos{27r{ki — k2)) cos {27r{ki + k2))\ dkidk2. 
2 Jo Jo 

If X = e^'^^^i and y = e^""^^^ , then 



1 rl 



Thex{z) =- 



log 



JO 



k- (x^ + x + + — 

^ y ^ J \ X y 



-2y X , y 

-2\ ( X , y 



dkidk2 



= -m ( k — ix^ + X ^) ( — h — ) ( X y H 

2 V \y xj \ X y 

We are now dealing with a Mahler measure. Let y — )• -. Then 



Thex{z) =^m [k - \x^ + 



1\ f X^ V 



2 ^' \ \ '' ' T.2 / V „ ' ^2 ' I 2/ + 



x^ J \ y x^ 



1 



m ((a;4 + l) (y^ + 1) (^^ + y^) _ kxS'^) . 



Finally let (x , y ) — ?• (x,y), and we arrive at 



= -m ((1 + x) (1 + y) (x + y) - /cxy) 
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Therefore, for — 1 < 2; < 1, we conclude that 



n,,{z) = ^g(^-l]. (27) 



When z = 1, we obtain 



Xhex = T/,ex(l) = ^5(8) = ((1 + x){l+ y){x + y)- 8xy) = 1^L_3(2). 

2 2 OTT 

The evaluation of g{8) is essentially due to Boyd, and a closely related calculation is 
described in [31, p. 31]. To recap Boyd's result, let us use (17) to obtain 

5 10 /5 4 1 1 

ff(8)=-log3--4F3(3,3^'i/;l 

= -m + 3xy) . 

6 

The Mahler measure n{k) := m (x^ + y^ + 1 — kxy^ also reduces to 4F3 functions [28]. If 
u) = g27r«/3^ then factoring the polynomial gives 

= -m [{x + y + l){x + uj'^y + u;){x + uy + u"^)) 
= -m(l + x + y), 

= ^L-3(2). 

The second step follows from additivity of Mahler measures, plus elementary changes of 
variable. The final step requires an evaluation due to Smyth [6]. 
Similarly, for the triangular lattice, 

TUz) = log 6 - - / ^^^^^ dz (28) 

7tJ z{3-z)y/{3-z){l + z) ^ ' 

where k = - — ^'^'^ ^'^'^ (^)' 



Xtri = log 6 







6V3K(A;) 1 



7rz2(3- z)V(3-z)(l + z) z 
3\/3 1111 ,1 

1.615329736 



(29) 



^^(5/6) \ 

52 + 72 112 + 132 + • • 27r j 



Now we would like to evaluate the function 



Ttri{z) = log 6 + j J ^°^{^~ ^ ^1 ^2 + cos(A;i + ^2))^ dkidk2 
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Suppose that < z < 1. Then we can insert an absolute value inside the logarithm, and 
we are again dealing with a Mahler measure. 



Ttri{z) =log6 + ^y j log 



(cos ki + COS k2 + cos(A:i + ^2)) 

Z O 



1 r-l 



JO 



log 



z 



+ 2 (cos(27rA;i) + cos(27rA;2) - cos(27r(A;i + ^2))) 



dkidk2 
dkidk2 



=m 



1 



6 1 1 

— hxH V y -\ xy 

z X y xy 



Now let {x,y) 1— )• [xy, |), to obtain 



1 , y , a; 



Ttri{z) =m [- + xy-\ \ \ y ^ 



=m 



=m 



xy X y 
+2+ ( x+ 



+ 2+ U + 







1 








X J 


y 






1 


^) 


[x^ 








X 



Finally, we make the substitution (x,?/) 1— )■ (xy, — f), and the formula becomes 



rf„(z) =m ( ^ + 2 - ( ^ + ^ ) ( X + 3 ) ( y + f 
=m ( (x + l)(y + l)(x + y) - ( ^ + 2 ) xy 



5 2 + 



In summary, we have the following result: 



Ttri{z) = g\2^- \ , 



6 



whenever z G (0, 1]. When 2; = 1 we obtain 



(30) 



\tri = Ttriil) = g{8) = 1^L_3(2), 

after appealing to evaluation of g{8) described above. In [16] Glasser and Wu give two 
other approaches to the evaluation of the spanning tree constant for the triangular lattice, 
including its expression in terms of the Clausen function, 



CH0) = Y^ 

n=l 



sin (n6 



It is Xtri = ^Ch (f ) . 
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From the connection between the LGFs, and the known expressions for their coefficients 
[18], we have the following identites: 

^-^ \ m I m \ ^ I TT 

m>0 



This is essentially due to Ramanujan. An equivalent version appears in [4] on page 40. It 
can also be readily proved by expressing both sides in terms of hyper geometric functions. 

An interesting equality between lattice sums arises from the honeycomb-triangular 
duality, Xtri/^ = Xhoney It follows that 

n>0 ^ ^ 3=0 ^ 
n>0 ^ ^ j=0 ^-^^ k=0 \ / \ / 

5.2 Three-dimensional lattices 

For the d-dimensional hyper-body-centred cubic lattice, the spanning tree generating func- 
tion is 

t!^^^{z) = log2'^ + ^ j\ki ■ ■ ■ j\kdlog (^^ - cos{ki) ■ cos{k2) ■ ■ ■ ■ cos{kd)^ (31) 

= log2'^H — -y / dki ■ ■ ■ / dfcfi log (1 — z cosf/ci) • cos(A;2) • • • • cosf/cd)) — log(z) 
Jo Jo 

= ''M2)-log(.)-i|^;(^)' 

2 

= dlog(2) - log(z) - ^d+2Fd+i [^'ll'-y-, 
So the spanning tree constant Ac is 

A^(d) = diog(2) - i ^ 1 ' (i^y = diog(2) - ^.+2F,+i 1 

a result first given by Chang and Shrock [12]. We give the results for d = 2, 3, 4 below, 
calculated to 50 significant digits almost instantaneously by Maple. 

A^^^(2) = 1.1662436161232751205535378258735796754562646159433 = iG/ir (32) 
A^^'=(3) = 1.9901914182719407717105190854333649929453600034709 
A^^'=(4) = 2.7329575354773621769814874419935610996620191115587. 

The result for the three-dimensional case also follows immediately from the Mahler mea- 
sure identification given by eqn. (19), so that 

A^^^(3)=p(64) = i/2(64). 
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Here G = 0.915965594177219 ... is Catalan's constant, and the result for d = 2 is well 
known as the square lattice spanning tree constant. For d = 3 and d = 4 the results lie 
outside the upper and lower bounds given by Felker and Lyons [13]. In both cases their 
upper bounds lie below the correct result. 

Because of the factorisation property peculiar to this lattice, the result is so simple 
that it has been derived without reference to Mahler measures. However it could have 
been derived that way, as the relevant Mahler measure is given by equation (18) above. 
The value of the spanning tree constant is obtained by evaluating the STGF at z = 1, and 
this evaluation is given by Samart in terms of L functions [36] as 

64 

A^-(3) = -L(r/(4r)6,3). 
Here r/(r) is the Dedekind eta function. 



5.3 Hyper simple cubic lattice 

The bcc lattice STGF above was the simplest to obtain. The corresponding s.c. result, 
as we show below, is the most difficult. For the ti-dimensional simple cubic lattice, the 
spanning tree generating function is 

Tf (z) = log(2d) + 1^ r dOi--- r dOdlog Q - i[cos(0i) + • • • + cos{ed)]\ (33) 



f'TT 

/ d9i-- 


■ £ dOdlogl^ 


Jo 








/ dOi-- 


■ £deaiog(^ 


JO 





oo 



1 1 I 

log(2d)-log(.)-ij:i(^) af^ 



1=1 

where the corresponding lattice Green function is 



1=0 

In particular, one has (see [17] for the 4-dimensional result): 

I /r\ 2 /oA 2 



(2) 



my- {'I 



j=0 



I 



,3,_/2!\„/A72A_/2!Y,^,A,-,,-, 



/ j=o \jy \J y \ 

It is a straightforward matter to sum the series (33), with z set to 1, to any reasonably 
desired accuracy with one's favourite computer algebra package, and observing that the 
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error made in stopping the summation after n terms decreases exponentially with n, we 
readily found 

Af (2) = T2{1) = 1.166243616123 . . . , 



A- 

Af 



(3) = r3(l) = 1.673389302970 . 

(4) = T4(l) = 1.999707644517. . 



Joyce and Zucker [25] have made a thorough study of this problem for the case -2 = 1, 
and have reported the value of related integrals for the case z = 1 to some 55 significant 
digits. The result in the two-dimensional case is exactly known, as given above, and we 
show below that we can evaluate the spanning tree generating function, and constant, 
exactly in the three-dimensional case also. In [42] Tzeng and Wu have also studied the 
problem of spanning trees on hybercubic lattices, as well as on non-orientable surfaces. 

From eqn. (18) we see that the simple cubic lattice STGF can be readily expressed as 
the logarithmic Mahler measure 

This is, in terms of a series expansion. 



Tsciz) = log(6) 



1 1 

2 ^ / 

1=1 



36 J 



E 



(34) 



Now we have from [32] 



oo 

E 

1=1 



u 



v9(l + 



u) 



^2V 



5(3 + 



u] 



/ 5 3 7 1 1 
-^Fa 4'2'4'1'1: 
\ 2,2,2,2 



'E 

j=0 

256^3 



27(1 + u)5 ^ 
(3 + n)3(l + 3u)^ 



(35) 



4 ' 2 ' 4 

"""" '9(3 + u)4 



+ 



4it 



15(1 + 3^)' 



l5 4 ^ 2,2,2,2 ' 9(1 + 3-^)4 



Unfortunately the argument of the second hypergeometric function is 1 when u = 1/9, so 
this expansion is only valid for u G [0, 1/9], whereas we require z G [0, 1] which corresponds 
to u G [0, 1]. 

That is to say, we would like to extend the formula (35) to all values of u G [0, 1]. The 
identity is valid when u lies in a neighborhood of zero, and it holds on the positive real 
axis if ti G [0, \]. It fails when u> because the argument of the second hypergeometric 
function intersects a branch cut on the interval [1, oo). It is easy to see that Q^^^g"^4 = 1 

when n = i. In general, we can analytically continue (35) along a ray from n = 0, until 

we reach a value for which either g^-^^g"-^4 G [l,oo), 

9^3+") 4 ^ [1)Cxd). When we cross a 
branch cut, it is necessary to correct an identity with additional terms, which are usually 
related to Meijer G-functions. 



16 



For numerical purposes, let us note that the following formula 



OO 

^ n V9(l + n) 

71=1 ^ 



2n 



n 



E 

A:=0 



2k 



n 



2, 

= 7 log 
5 



27(1 + 



(3 + n)3(l + 3u)7 5(3 + ^)4 

2 



+ 



,1,1. 



256^3 



2,2,2,2 '9(3 + ^)4 



(36) 



1 /■"' \2Fi ( 4j4;t 



10 



1 



t{l-t) 



-(1 - 2t)dt, 



holds if 



9r 



(l-r)(8 + r)' 



r(8 + r)" 



+ 20r -r2)^' 



for ah r G [0, 6\/2 - 8]. This restriction on r implies that u £ [0, 1] and w £ [0, ^±|^] w 
[0,0.97]. The integral on the right is analytic when \w\ < 1, so it is easy to see that 
formula (36) is valid on a somewhat larger domain than formula (35). Formula (36) is 
an intermediate step in the derivation of (35), but is not actually stated in [32]. The 
two identities coincide if w £ [0, ^], because Clausen's identity allows us to perform the 
integration [9]: 



w 2 



t{i-t) 



'1 - 2t)dt 



w 3J 



F2[™-At{l-t))-l 

t{l-t) 



(1 - 2t)dt 



4u.(l-t«) 3F2 4'2'4;tM - 1 

^ ^ dt' 



^'^wil - w)^Fi (^^^I;4w;(l - w) 



Since Aw{l —w) = g^f^^yri we recover (35). Clausen's identity fails when w > 1/2, so we 
need an additional method to simplify (36). 

We simplify (36) by using modular parameterizations for hyper geometric functions. 
Consider the following g-series: 



G{q) := K 



log((?) + 240 5^n2log(l-g'^) 



n=l 



and note that 



G{q)+Gi- 



9G(<?2) -4G(g^). 



(37) 



(38) 



Formula (38) is proved in [32, Thm. 2.3]. Recall that /2(ii) and /^{u) are defined in eqn 
(20), and define the sum 



log('u) - 



(1/u)^" f2n 



n=l 



2n 



n 



2k\ In 



fc=0 



(39) 
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Results of Bertin [2] and Rogers [32] show that if l^l is sufficiently small, then 

aiihiq)) = - ^G{q) + ^G{q') - ^G{q') + ^G{q'), (40) 

f2{s2{q)) = - j^G{q) + ^G{q% (41) 

h{sM)=-lG{q) + lGiq''), (42) 



where 



^^^^^ ^ / r]{q)ri{q^) \ ^ ^ ( v{q)viq^) 



S2{q) 



\viq^)r]{q^) J \viq'^)viq^ 

( vHq') 
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KviqMq"^) 



and 



oo oo 



ri{q)=q'"'\{{l-qn= E (-1) V^^+^^'Z^^ 

n=l n=— oo 



We note that Rogers gave a slightly different version of (41), and that this form of the 
identity was first observed by Samart [36]. Rodriguez-Villegas and Stienstra proved similar 
modular expansions for two- variable Mahler measures [39], [31]. Using formulas (38), (40), 
(41), and (42), we deduce the following five term relation: 

9l{tl{q)) = - \f2 {S2 i-q)) + (S4 {q^)) + {sa {q')) 

4 20 20 (43) 

+ Y^h{s.{-q'))-\h{s.{q')), 

which holds for qm. a. neighborhood of zero. There are many additional relations between 
the hypergeometric functions, but this is the simplest formula we were able to find which 
holds for g G fo, e^"""^] . We require the identity to be valid on a region which includes 

q = e^'^V^^, because the s.c. constant is given by (71(6) = gi (^i ^e~'^V^^^^. Formula 

(43) holds in a neighborhood of g = 0, and it can be extended to the desired region, 
because the arguments of the hypergeometric functions never intersect the line [l,oo\ 
For instance, we can use a computer to check that g2(^q^ [1? c«) for all q € 0, e^'^^f^ . 

Equation (43) simplifies if we observe that the functions {ii(g), S2{—q)^ sa{()^\ • • • } are 
algebraically dependent. To obtain explicit relations, we need modular equations from 
Ramanujan's notebooks [4]. Suppose that q and /3 are given by 

a := 16 , „. — , p := 16- 



^24(^2) ' ^ . ^ rf^{qi3) ■ 

It is known that a and /3 are algebraically related if j is a positive integer, and we call such 
relations jth degree modular equations (this is equivalent to Berndt's definition [5, p. 212]). 
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When j = 3 we appeal to [5, p. 230], to see that 

(a/3)i/4 + {(l-a)(l-/3)}i/4 = l. 

Ramanujan parameterized the third-degree modular equation in terms of rational functions 
[5, p. 230]: 



We easily solve (44) for p: 



p — z- 



By the results in [5, pg. 124-126], we also have 



\v{q) J 16(1 -a)2' Kviq'^) J 256(1 - a) ' 

v{q')V'_ a(l-a) f v{q') _ {l - VT^f VT 



vi~q) J 16 ' \'ni~q'^) J 



a 



(45) 
(46) 

, _ (47) 

Mq^)J 1024a2Vr^ 

If g I— >■ in either (45), (46), or (47), then the identities remain valid as long as we replace 
a with /?. Applying (44), (45), (46), and (47) leads to: 



s2{-q) 

i-q') 



1-/3/ \ct J \1 — a J \/3 
2(1 + p + p2)(i + 4p + p2^ 

"(l-p2)^p(2 + p)(l + 2p) 
16(1 - a)2 



a 

4 



16(2 - a 
a2(l — a) 
16(1 + /3)^ 
^/3(l-/3)2' 

4 (1 - ^/T^)'^ (2 - a - QVT^Y 

4 (1 + (2 - a + eVT^^"^ 



After substituting the above parameterizations into (43), we can drop the dependence 
on q, and simply regard a and /3 as functions of p. It is possible to prove that p = 
i (2 + 3V2 - 2V3 - Vg) when q = e"^^/^ . The derivation is an exercise in manipulating 
singular moduli [7, p. 183], and is easy to check numerically. 
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Theorem 1. The following identity is valid: 

'2(1 + p + + 4p + p2)\ 



91 



(l-p2)^p(2 + p)(l + 2p) 



4 \ a J 20 V«^(l-a)y 20 V/3(l-/3) 



1 / 4(1- vr^)'(2-a-6^/r^)'\ ^^^^ 



10'" \^ a4^/i^ 

1 /4(i + vr^)^(2-a + 6^/r^) 

7/' 



5'" V a'^VT 



a 



provided that 



l + 2pj ' " Vl + 2p 

andpG [0, 1(2 + 3^/2-2^/3-^/6)] « [0, .16452 ... ]. 

If we wish to compare formulas (36) and (48) numerically, it is sufficient to set r = 
in (36). The left-hand side of (48) equals 51(6) when p = i (2 + 3\/2 - 2^/3 - ^6) . 
Since gi(6) = Xsc, we obtain 

^ - ^fj 16(1 -«)^ \ , 9^^/ 16(2-a)n ^ 3 16(1 + /?)^ 



4'^V « J 20' Va^(i-") y 20' V/?(i-/3) 

1 / 4(1- Vr^)^(2-a-6\/r^)^\ 



10' \^ a^^/T^ J 

1 /4(l + vT^^^(2-Q + 6vT^)^\ 
- -/4 — ^ \ \ ^— = 1.673389302970196732283. . . . 

(49) 

We defined f2iu) and /4(n) in (20) and (20), and a and /3 equal 

a =35 — 24\/2 - 20\/3 + 14\/6, 
/3 =35 + 24V2 - 20V3 - UV6. 

Since all five of the hypergeometric functions are real-valued, we can drop the "real part" 
notation from the hypergeometric functions in (20) and (20). Note that a slight simplifi- 
cation arises if we observe that the argument of the third term above simplifies to 2304. 
No other terms obviously simplify. 

5.4 L-function formula for the spanning tree constant. 

It is possible represent A^c in terms of L- functions of modular forms [34]. We have 
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where q = e^'^*^, and 



2 

n,fc=— oo 

We can also represent /(r) in terms of eta functions: 

^ _ 7?(r)2r? (4t)S (6r)3 r?(r)4r? (4rf (6t)S (8rf 
r?(2r)%(8r)2 77 (2r)% (3r)% (12r)2 

^ r?(2r)^?7(8r)%(12rf 7?(r)2r? (2r)=^ (8r)^ (12r)^ 



r?(r)2r/ (4r)^ r? (4r)% (6r)^ r/ (24r)^ 

Bertin [3] has obtained some similar results along these lines. 

5.5 Diamond lattice 

For the diamond lattice, using the abbreviated notation q = cos(A;i) i = 1,2,3, 

Tdiam{z) = log 4 + ^ j j j ~ ~^ ^^'^^ "^^^^ ~^ ^203)^ 

- .o.-iE(rkc)^a)(::r 

n>l i=0 ^"'^ \ J / \ J 

The argument of the logarithm above is non-negative for £ [0, 1], so can be replaced 
by the logarithm of the modulus for z in that range. From eqn. (18) one then immediately 
recognizes the integral above as the Mahler measure f{16/z^) — 41og(2), so that 

Tdiam{z) = f {16/ z^)-2log{2). 

Substitution into eqns (19) and (20) and straightforward algebraic manipulation then gives 

Td^an^{z) = ^ log(2) + ^ log(4 - Z^) _ 1 log(l - ^2) 

5 5 ID 

^2 „ A 3 4 , T -27z^ 



40(1-^2)3= " 2,2,2,2 '4(1-^2)3 

4z^ ^ A 3 4 1 1 27z^ 



5(4-^2)3° ^ \^ 2,2,2,2 ' (4_^2)3y " 

This expression is well defined for \z\ < 1. As z — t- 1, both a logarithmic term and the 
first hypergeometric term are singular. We wish to calculate 

1 , / 9x ^ /s 3 4 , 1 -272;2 \ 

hm -- log(l -z)- ^^^^—^,.F, [^^^^ ■ ^^-^^ j . (51) 
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If we express the 5-F4 function as an integral, then the hmit becomes 



112 



lim— — log — h — / ^ an + — / ^ —au 



1 30 4 30 7o ^ 30 

1 27 1 /■! 3^2 - 1 1 /-oo 3F2 ( S'f'3;-n' 

= log \ / — du H / — <-du 

30 ^ 4 30 io u 30 7i u 

2 , 
=— log2. 

15 ^ 

Mathematica can evaluate both definite integrals in terms of 4-F3 functions. When the 
expressions are combined, all of the 4-F3 functions cancel out, and we are left with a 
multiple of log 2. 

Thus we obtain for the diamond lattice spanning tree constant 



\d^a^ = log(2) + - log(3) - — 5F4 ( 3'5.3.^^.^; 1 ) = 1.2064599496517629 . . . 



5 3 4 



Prom eqn (2.25) in [32] we have the elegant result giving the spanning tree constant 
in terms of an L-function, 



24 

\diam = — L(??(2r)=^7?(6r)3,3) 



where r/(r) is the Dedekind eta function. 

For the 4d diamond lattice, since the LGF is identical to that of the 4d hyper-cubic 
lattice, and the two lattices have the same co-ordination numbers, it follows that the 
spanning tree constant of the 4d diamond lattice is the same as that of the 4d hypercubic 
lattice, given above. 



5.6 Face-centred cubic lattice 

After the completion of this work, we became aware of the work of G S Joyce [23] who 
studied both the STGF and LGF on the f.c.c. lattice making use of the relevant Mahler 
measure given by one of us in [32]. Accordingly, we have shortened this section and just 
given brief details of our derivation. A more extensive and thorough treatment of the f.c.c. 
case can be found in [23]. 

The f.c.c. spanning tree constant was first given by Shrock and Wu [37], but there 
was an error in the structure function that they used. This was corrected by Chang and 
Shrock [12] who correctly gave \fcc ~ 2.41292. The spanning tree generating function is 

I 1 

Tfcc{z) = log 12+- — -:7 / dki--- / dhi\og{l/z — (cos fci cos A;2+cos A;2 cos fcs+cos fca cos 
(27r)'^ J__^ 3 

Prom eqn (18) we immediately identify Tj^c as /(4 + 12/z). Substitution into eqns (19) 
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and (20) and similar algebraic manipulation as in the diamond case above then gives 

Tf,,{z) = 21og(2) + ^log(3)-log(z)-^log(3 + z)-ilog(l-z) (52) 

5-^4 I 3^2^ 3^ 



135(z-l)3^ ^ V ^'^'^'^ ' (^-1)^ 
2z\z + 3) Aj_4 i i_z2(2 + 3) 



135 M 2.2,2,2 ' 4 



As above, taking the limit as z — t- 1 we obtain for the fee lattice spanning tree constant 
A/cc = 2Xdiam = ^ log(2) + ^ log(3) - r^^Fi (t^ll'J-^ l) = 2.4129198993035259 . . . 

As for the diamond lattice case, from eqn (2.25) in [32] we have the result in terms of 
L- functions, 

48 

A/cc = ^^(^(2T)^?(6r)^3). 

Finally, for the four-dimensional fee lattice we find by integration of the relevant LGF 
[17], A/cc(4) « 3.14567.... 



6 Integral identities and hyper geometric identities 

Having shown the connection between STGFs and LGFs, we can immediately write down 
a number of integral identites which follow directly from the fact that the integral of the 
LGF is the STGF. As two- and three-dimensional LGFs can be expressed as complete 
elliptic integrals of the first kind, and as the square of complete elliptic integrals of the 
first kind respectively, while STGFs are expressible in terms of hypergeometric functions, 
it follows that we can express the integrals of complete elliptic integrals and their square 
in terms of higher order hypergeometric functions, while differentiating the STGF result 
relates higher order hypergeometric functions to 2F1 hypergeometric functions, or their 
square. 

In this way a large number of identities can be produced. We have not yet determined 
how new and/or useful these results are, so we just give a couple of examples at this stage, 
thereby making the procedure clear. We will investigate the full family in greater detail 
subsequently. A simple, previously known, result is given above in equation (24). Less 
well-known, or possibly unknown results follow in the case of three-dimensional lattices. 
The simplest situation arises with the body-centred cubic case, by combining eqn (8) and 
the equation above for A|!'^'^(3), giving the integral identity 



MfK(M]^-l _ 1 



"^^^ i6'^HT2" ,2'^ 



1,1, 



3 



where 

2 2 



fc2 = i_iyr3^ 
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More generally, 



\ dz = \ogz + -5F4 (''2£':;2^; 



Similar identities can be obtained for all the other two- and three-dimensional lattices. 

For example, by combining eqn (8) and the equation above for A^'^(3), we now have 
the integral identity 



1 / 4(1- VT^)^(2-a-6^/r^)^\ 



(53) 



/4 



1 /4(l + \/r^)^ (2-a + 



where 



with 



^^^^ 



(l_^)3(l + 3^): 



e = (1 + ^/l-Z^)~^l\l - 71^12/^)1/2 

and several other identitites also follow immediately. 

As another example, by combining eqn (8) and the equations above for A/cc, we now 
have the (possibly new) result: 

r ii'mlL) [fK(fc3)]'-l 48 , , 

Jo Z 

where 

3 (l_^)3(l + 3^)' 

and 



e = (1 + + VTT3^). 

An alternative expression for A/cc is 

v^^-o..-Ei(i^)"t(;)M.-E(0\l^)(^-f). 

n>0 ^ ^ j=0 ^-^^ fc=0 \ / \ / \ / 

Note too that A^ = Xfcc/'2, giving rise to the identity 

EK.^)E„(;).-.-E(Ott)(-f)-.o.! 
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Going in the opposite direction, for the square lattice one can, for example, differentiate 
equation (23), giving 

-^-T = 1 + ^4^3 %^V;^ +T7T^4i^3 %%^o^;z' (54) 



dz 4 " " V 2,2,2 ' ; 128" " V 3,3,3 

= i + T + ^^^H3^3^'")=^^HT'" 

The second equality follows by algebraic manipulation of the sum of the two higher order 
hyper geometric functions, and the final equality follows as the r.h.s is the square lattice 
LGF. 

A corresponding calculation based on the results above for the body-centred cubic 
lattice STGF and LGF gives 



Again, similar identities can be obtained for all the other two- and three-dimensional 
lattices. 



7 Connection between spanning trees and the Ising model 
and dimer coverings. 

There is a close connection between the spanning tree constant A and the free-energy of 
the Ising model at the critical temperature. Unfortunately, this seems to be true only 
for planar lattices (otherwise we would have some clue as to the structure of the free- 
energy of the 3-dimensional Ising model) . The Onsager solution for the free-energy of the 
square-lattice Ising model is 

F{v) = log (^ ^ ^^2 ) + 8^ / / ^og[{l + v'^)'^ -2v{l-v'^){cosx + cos y)]dx.dy (56) 

where v = tanh(jj^). At the critical temperature, v = Vc = — 1, this simplifies to 

F{vc) = ^log{2) + ^ J J log[4-2(cosx + cos2/)]dx.dy = (A,^ + log2)/2. (57) 

If bn denotes the number of distinct dimer coverings of an n x n square lattice (with 
n even), then, as shown by both Kasteleyn [26] and by Fisher and Temperley [41], 

n/2 n/2 . 

6» = 2"/^nn(<=»^;rjT + ~=' 



j=l k=l 

In the infinite lattice limit, 



n + 1 n + 1 



1 1 /""' G 

lim — log6„ = 7j / / log[4 + 2(cos X -|- cos y)](ix.(i?/ = — , 

n-5>oo n lovr^ 1^ 1^ vr 
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where G is Catalan's constant. This is very similar to the expression for A^g, apart from 
a (critical) sign change. However Temperley [40] pointed out that if one considers dimers 
on a [2n — 1) x [2n — 1) site lattice, with one boundary site removed, the relevant sign 
changes, and the integral agrees with that for the spanning tree constant. Later, Tzeng 
and Wu [43] evaluated the dimer generating function for the one vacancy case, and showed 
that it was independent of the location of the vacant site. The Temperley bijection does 
not apply to cylindrical lattices however. For that situation, the relevant dimer generating 
function has been evaluated by Wu, Tzeng and Izmailian using Pfaffians [46]. 

In [39], Stienstra points out a connection between the partition function of certain 
dimer models and the L-functions of their spectral curves. This follows from his obser- 
vation that the partition function per fundamental domain of a dimer model, as given 
by Kenyon, Okounkov and Sheffield [27] is in fact the logarithmic Mahler measure of the 
characteristic polynomial of that dimer model. Stienstra demonstrates the connection 
explicitly for three different dimer models. 

8 Conclusion 

We have derived the spanning tree constants for all the usual three-dimensional lattices 
in terms of 5F4 hypergeometric functions. This has been made possible by establishing 
the connection between the integral representation of the spanning tree constant and the 
Mahler measure of an approporiate Laurent polynomial, which is closely related to the 
structure function of the underlying lattice. We have introduced the notion of a spanning 
tree generating function, which gives the spanning tree constant as the value of the STGF 
at z = 1. We show the simple connection between the STGF and the lattice Green function, 
and so by comparing known results for the LGF with known and new results for the STGF, 
we are able to derive a number of integral and hypergeometric identities. We have also 
expressed all of the spanning tree constants for both two- and three-dimensional lattices 
in terms of Dirichlet L-series. 
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